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In this study, we extend the results of previous combined numerical and experimental investigations of 
an axisymmetric laminar diffusion flame in which difference Raman spectroscopy, laser-induced fluores- 
cence (LIF), and a multidimensional flame model were used to generate profiles of the temperature and 
major and minor species. A procedure is outlined by which the number densities of ground-state CH 
(X 2 fT), excited-state CH [A 1 A, denoted CH*), and excited-state OH (A 2 £, denoted OH*) are measured 
and modeled. CH* and OH* number densities are deconvolved from line-of-sight flame-emission mea- 
surements. Ground-state CH is measured using linear LIF. The computations are done with GRI Mech 
2.11 as well as an alternate hydrocarbon mechanism. In both cases, additional reactions for the production 
and consumption of CH* and OH* are added from recent kinetic studies. Collisional quenching and 
spontaneous emission are responsible for the de-excitation of the excited-state radicals. 

As with our previous investigations, GRI Mech 2.11 continues to produce very good agreement with 
the overall flame length observed in the experiments, while significantly under predicting the flame lift- 
off height. The alternate kinetic scheme is much more accurate in predicting lift-off height but overpredicts 
the overall flame length. Ground-state CH profiles predicted with GRI Mech 2.11 are in excellent agree- 
ment with the corresponding measurements, regarding both spatial distribution and absolute concentration 
(measured at 4 ppm) of the CH radical. Calculations of the excited-state species show reasonable agree- 
ment with the measurements as far as spatial distribution and overall characteristics are concerned. For 
OH*, the measured peak mole fraction, 1.3 X 10 compared well with computed peaks, while the 
measured peak level for CH*, 2 X 10‘ 9 , was severely underpredicted by both kinetic schemes, indicating 
that the formation and destruction kinetics associated with excited-state species in flames require further 
research. 


Introduction 

CH has long been recognized as a key reactant in 
NO* formation through the prompt NO mechanism. 
Given that CH is a short-lived trace species that ex- 
ists in a narrow spatial and temperature region 
within a flame, its concentration and spatial distri- 
bution are very sensitive tests of the detailed chem- 
ical kinetics needed to model pollutant formation. 
The oxidation of CH plays a central role in the pro- 
duction of chemically excited OH (A 2 Z, denoted 
OH*), which emits in the ultraviolet. This ultraviolet 
emission has been suggested as a measure of the 
final steps of the CH* reduction chain [1]. Chemi- 
cally excited CH (A 2 A, denoted CH*) is responsible 
for the blue light in low-soot flames and may provide 
insight into the C 2 reaction chain [1]. Despite the 
prevalence of CH* and OH* chemiluminescence, 
little quantitative work has been done either exper- 
imentally or computationally in predicting the ab- 
solute concentrations of these species. In this study, 
quantitative measurements of CH, CH*, and OH* 
are performed and comparisons are made with com- 
putational predictions. 


The flame under investigation is a lifted axisym- 
metric laminar diffusion flame, which has been char- 
acterized previously both experimentally and com- 
putationally [2-5]. The fuel is nitrogen-diluted 
methane surrounded by an air coflow. Experimen- 
tally, temperature and major species (CH 4 , N 2 , 0 2 , 
H 2 0, CO, C0 2 , H 2 ) concentrations were measured 
simultaneously with Rayleigh and Raman scattering 
[2]. Laser-induced fluorescence (LIF) measure- 
ments were performed to measure number densities 
of minor species. Quantitative, linear LIF measure- 
ments were made for OH [3] and NO [4], and qual- 
itative measurements of CH have been made [3]. 
Modeling work has employed different kinetics 
schemes, including a 26-species C 2 hydrocarbon 
mechanism [3] and GRI Mech 2.11. Both produced 
excellent agreement for temperature and major spe- 
cies [5]. Computed peak concentrations for NO and 
OH were within 30% and 15%, respectively, of their 
measured values. 

In the following sections, the experimental config- 
uration is described, and the details of the measure- 
ments of CH, CH*, and OH* are presented. The 
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computational model and the various kinetics 
schemes used to predict the measured species are 
then described. Finally, the experimental results are 
compared with computations based on different ki- 
netics schemes. 


Burner Configuration 

The burner used in this experiment consisted of a 
central fuel jet (4 mm diameter) surrounded by cof- 
lowing air (50 mm diameter). The fuel was com- 
posed of 65% methane diluted with 35% nitrogen 
by volume to reduce soot, and the plug flow exit 
velocity of both fuel and coflow was 35 cm/s. This 
produces a blue flame roughly 3 cm in length with 
a lift-off height of 5.5 mm. Complete burner speci- 
fications are given elsewhere [2]. The burner was 
mounted on a stepper motor to allow measurements 
to be taken at different heights. 


Laser-Induced Fluorescence Measurement 
of CH 

The third harmonic of a Nd:YAG laser, operating 
at 10 Hz, pumped a dye laser containing Coumarin 
440 dye. The R(7) line in the (0,0) band of the A-X 
system, near 426.8 nm, was selected for excitation. 
The dye beam was split for power measurement, at- 
tenuated, shaped into a sheet, and passed across the 
jet centerline. Typical energies were 1 /uj per pulse. 
The laser sheet dimensions were measured to be 5.5 
mm X 300 //m. The measured line width Av of the 
beam was 0.16 cnr 1 , and the pulse duration was 10 
ns, which corresponds to a spectral density of order 
10 4 W/(cm 2 cm -1 ), well below the saturation value 
[6]. Imaging was done with a cooled charged-cou- 
pled device (CCD) camera and a lens-coupled image 
intensifier. The CH fluorescence was isolated with 
an interference filter that transmitted from 400 to 
440 nm. An f/2 camera lens collected the CH fluo- 
rescence, and a pair of f/1.4 camera lenses focused 
the light from the back end of the intensifier onto 
the CCD chip. The imaged pixel volume was 30 X 
30 X 300 //in 3 . The laser was set to 426.777 nm to 
record on-resonance images (Z on ) and to 426.671 nm 
for off-resonance images (/ 0 n)- A flame luminosity 
image was also taken with the laser turned off (J| nin ). 
Typical images were integrated over 6000 laser 
shots. The on-resonance image contained LIF and 
Rayleigh scattering on top of the flame luminosity 
background, while the off-resonance image had only 
Rayleigh scattering and flame luminosity. Because 
the laser energy differed by a small known amount 
between the on- and off- resonance measurements 
(E on and E ( , B ), the final LIF images were created as 
follows. Sjjj? / on hum (h>H hum) 
Measurements were made at each of eight heights 


above the burner surface, and these images were 
later tiled together. 

To convert the measured LIF signal into a quan- 
titative concentration measurement, a number of 
calibrations and corrections must be made. For a 
two-level model, the LIF signal per pulse in the lin- 
ear regime is given by 
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where b\ 2 is the absorption rate, Euf is the laser 
energy per pulse, Av is the laser line width, A up is 
the cross-sectional area of the laser beam, n (:11 is the 
number density of ground-state CH ,/ B is the frac- 
tion of the ground state in the state being pumped, 
Eis a dimensionless overlap integral, A 2 j is the spon- 
taneous emission rate, @21 * s collisional quench- 
ing rate, Vj IK is the LIF pixel volume, Q is the solid 
angle over which light is collected, £ is the efficiency 
of the collection optics, and ?/ is the detector effi- 
ciency in counts per photon [7]. The quantity that 
we are interested in determining is Uc:m- The factors 
in equation 1 must be measured directly or deter- 
mined from the literature. The quantity A 2! /( A? ] + 
@ 2 i ) represents the fraction of excited molecules 
that emit a photon and is called the fluorescence 
yield, 0. The total collisional quenching rate for the 
excited state is @ 2 1 — Scjk,, where c, is the concen- 
tration of species i and k i is the collisional quenching 
rate for species i. Values for k { (which are tempera- 
ture dependent) w’ere computed from the functional 
form given by Tamura et al. [8]. Concentration pro- 
files of the major species (CH 4 , N 2 , 0 2> H 2 0, CO, 
C0 2 , H 2 ) were computed previously [5] and used 
for these calculations. The temperature in the thin 
region where CH is present is nearly constant at 
1900 K. This quenching calculation resulted in a spa- 
tially constant fluorescence yield 0 of 1/176. In this 
flame, N 2 was responsible for more than 60% of the 
quenching. Unfortunately, the CH/N 2 quenching 
rate has not been measured above 1300 K making 
the extrapolation to flame temperatures a source of 
uncertainty. The Boltzmann factor, / B , is calculated 
based on the molecule probed, temperature, and the 
excitation/detection scheme used. Given our choice 
of the R(7) line,/ B = 0.07. Spectroscopic rotational 
and vibrational constants were taken from Huber 
and Herzberg [9], Eis a spectral integral of the over- 
lap between the CH absorption profile and the laser 
beam. The laser profile was fit to a Gaussian distri- 
bution centered at 426.777 nm with a full width at 
half-maximum (FWHM) of Av ~ 0.16 cm -1 , mea- 
sured previously. The absorption spectrum of CH (at 
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Fro. 1. CH mole fraction profiles determined In measurement, GUI Mech 2.11, and an alternate hydrocarbon 
mechanism. 


atmospheric pressure with resolution Ay) was cal- 
culated with the LIFBASE program [10] and re- 
sulted in r = 0.17. The calculated absorption profile 
was needed because aCH excitation spectrum could 
not be measured practically in the linear regime 
given the long integration times needed to obtain 
reasonable signal levels. A saturated CH excitation 
scan was performed, however, and compared well to 
a computed absorption profile for a larger line width. 
When evaluating the integral, both the computed 
absorption profile (with line width Av) and fit beam 
profile had narrow features. Small variations in the 
peak of the Gaussian, hence small changes in the 
laser wavelength, affect r. Given that the resetta- 
bility of the dye laser is ± 0.02 A, the uncertainty in 
r is estimated at 20%. 

Rayleigh Calibration 

Calibration of CH LIF is somewhat problematic. 
CH is highly reactive, preventing reference to a 
known concentration, and is present in a very thin 
(—500 fi m wide) region, making absorption calibra- 
tion difficult. However, Luque and Crosley [7] have 
demonstrated that Rayleigh scattering can be used 
on the same optical setup to relate the measured 
signal to an absolute light level and thus solve for 
the overall calibration product (Qerj) in equation 1, 
For scattering from a homogeneous gas, the Ray- 
leigh signal per pulse is given by 

S K = —j Qen (2) 

A R /iv 


where N is the number density of the gas flow used 
for calibration, E n is the laser energy per pulse pro- 
ducing the Rayleigh signal, V R is the Rayleigh vol- 
ume, (da/dQ) is the Rayleigh cross section, hv is the 
photon energy, and A H is the cross-sectional area of 
the laser beam. After Rayleigh measurements are 
made, the calibration product can be eliminated, re- 
sulting in a single expression for number density 
Given that the fluorescence beam occupies the same 
spatial region as the Rayleigh beam (Ar = A LIF and 
Vr = V UF ), we now have 

4*S UF JriVE R (^) 

" c;H = b, 2 E uv rf n 0hvS H (3) 

as an expression for absolute number density of 
ground -state CH. 

Rayleigh calibration was done on a flow of clean 
air with the same optical setup and beam dimensions 
as used in LIF. Although calibration can be done 
with a single Rayleigh image, measurements were 
made over a range of energies to verify the linear 
relationship between the signal and laser energy and 
to confirm a zero intercept. The measured calibra- 
tion was linear and had an intercept within a few 
percent of the origin. Energy per pulse was varied 
between 0.1 and 1.2 mj, with an integration of 300 
to 1500 shots. 

These calibrations and calculations, in conjunction 
with previous temperature measurements, resulted 
in a two-dimensional profile of CH mole fraction 
that can be seen in Fig. 1. The measured peak mole 
fraction was 4 ppm. 
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Chemiluminescence of CH* and OH* 

Measurements of chemically excited Hanie radi- 
cals, such as Cl I* and OH*, arc relatively easy to 
make hut have calibration difficulties similar to those 
of ground-state CH. Fortunately, the same Rayleigh 
calibrations can he applied to determine absolute in- 
tensity levels, and quenching corrections can he 
made to quantify the measured signals. The space 
and power constraints imposed by micro gravity ex- 
periments have generated a renewed interest in 
making such emission measurements as quantitative 
as possible. Although flame emissions were exam- 
ined spectrally decades ago and CH* has been rec- 
ognized as the primary source of the ubiquitous 
“blue light” in low-soot flames [11], the reactions 
that produce CH* are not well understood. Chemi- 
luminescent OH (OH*) is another common flame 
emitter whose kinetics deserve inv estigation. A goal 
of this study was to measure the absolute number 
densities of these flame radicals to assess the current 
state of CH* and OH* kinetics. 

OH* chemiluminescence peaks at 307.8 nm. Ex- 
cited-state OH measurements were made with a 
cooled CCD camera using an f/4.5 UY camera lens; 
the camera/lens system was placed 50 cm away to 
ensure a wide depth of field. A narrow bandpass UV 
filter was used (center 307 nm, 10 nm bandwidth). 
High emission signal levels were collected with 10 s 
integration times. Rayleigh calibration was done 
with a Nd:YAG-pumped dye laser utilizing sulfar- 
hodamine 640 dye. The 612 nm output was fre- 
quency doubled to perform calibration measure- 
ments on clean air. Beam energy in the UV was 
varied from 0.3 to 2.0 mj per pulse and integration 
times were 600 shots. 

For CH* flame emission, the A~A — > X 2 J7 tran- 
sition at 431.2 nm was imaged with the camera and 
lens system detailed above. An interference filter 
(center 431 nm, 10 nm bandwidth) was used to iso- 
late the CH* emission. Again, good signal levels 
were collected over 10 s integration times. The cal- 
ibration measurements were again performed with 
a Nd:YAG-pumped dye laser. During Rayleigh cali- 
bration, the energy of the 431.5-nm beam was v aried 
between 0.5 and 2.5 mj, with signal collected over 
600 shots. 

The flame emission signal, in detector counts, is 
given by 

s..„, = T A. 2i TV„ n N v „fiaj (4) 
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where t is the integration time, Y t . m is the pixel vol- 
ume of the emission signal, and A 7 L . m is the number 
of molecules that emit a photon. The spontaneous 
emission rate and integration time are known, and 
the emission volume can be determined directly. As 
with the fluorescence measurement, the calibration 
constants of the optical system, Qetj y are accounted 


for with Rayleigh scattering. Furthermore, a 
quenching calculation is needed to relate the num- 
ber density of emitting molecules to the total popu- 
lation of the chemically excited radical. To obtain an 
absolute number density of a chemically excited rad- 
ical. A 7 *, we recognize that A 7 * = N cu /0. Equations 
(2) and (4) can he combined with this to yield 

** 

N* = — (5) 

S^zV^hv 

where € is the length along the beam in the Rayleigh 
volume under consideration. 

Emission measurements are integrated through 
the collection optics along the line of sight. Appro- 
priate background images, taken for both CH* and 
Oil* with the flame extinguished, are subtracted 
from the raw emission signal. Given that our flame 
is axisvm metric and that the imaging optics are con- 
figured so that the magnification changes by only 17c 
over the flame width, we can recover a tw o-dimen- 
sional, in-plane intensity distribution proportional to 
number density with the use of an algorithm that is 
equivalent to a hvo-point Abel deconvolution [12]. 
After inversion, pixel volumes were determined to 
be cubes of side length 68 //m for both CH* and 
OH*. The quenching correction for CH* was per- 
formed using the cross sections and temperature de- 
termined for ground-state CH. This resulted in a 
peak CH* mole fraction of 2 X 10“ as seen in Fig. 
2. For OH*, quenching rates for OH were taken 
from Tamil ra et al. [8] and combined with major 
species and temperature computations performed 
previously [5]. OH* appears at a nearly constant 
temperature of 1900 K and is highly localized. The 
quenching calculation resulted in a fluorescence 
yield 0 of 1/327 that did not vary spatially. Water, 
the only important individual collider, is responsible 
for about 2/3 of the OH* quenching in this flame. 
Cross sections for OH have been measured over a 
large range of temperatures and have agreed with 
low-pressure quenching measurements to within 5% 
[8]. The measured OH* profile is displayed in Fig. 
3; the peak mole fraction was 1.3 X 10 ~ s . 

Computational Model 

The computational model used to compute the 
temperature field, velocities, and species concentra- 
tions solves the full set of elliptic two-dimensional 
governing equations for mass, momentum, species, 
and energy conservation on a tw o-dimensional mesh 
[13]. The resulting nonlinear equations are then 
solved on an IBM RS/6000 Model 590 computer by 
a combination of time integration and Newtons 
method. The chemical mechanisms employed were 
GRI Mech 2.11 [14] and an alternate hydrocarbon 
mechanism (used, for example, in Ref. [3]). 
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Fig. 2. CH* molt 1 fraction profiles 
determined by measurement, GRI 
Mech 2.11. and an alternate hydro- 
carbon mechanism. The ground- 
state CH profile is shown to illustrate 
the thin features of both the mea- 
sured and computed CH* profiles. 


Accurate computations of the CH, CH*, and OH* 
radicals pose a difficult numerical problem. Forex- 
ample, the concentration of CH can change by an 
order of magnitude within 0.1 mm. This requires 
that the adaptive grid be refined to an extremely 
small mesh size in the v icinity of high spatial activity. 


Kinetics Modeling 

To model chemiluminescence, the species CH* 
and OH* must he added to the kinetic mechanisms, 
along with a set of formation and destruction reac- 
tions with appropriate rate constants. All CH*/OH*- 
related rate constants used in this study are detailed 
in Table 1. 

CH* is produced chemically via the reaction of the 
ethynyl radical with monatomic and diatomic oxy- 
gei 


I) C 2 H + 0 2 CII* + C0 2 

II) C 2 H + O CH* + CO 

The rate constants have been reported as k } = 3.6 
X 10" 14 cm 3 molecule - 1 s" 1 and k n — 1.8 X 10 -11 
enr 1 molecule -1 s -1 [15]. The uncertainties associ- 
ated with these reaction rates are about 40%. De- 
struction reactions occur by spontaneous emission 
(CH* — > CH 4- h v) and collisional quenching. 
Quenching was modeled with seven different reac- 
tions, each involving CH* and a major chemical spe- 
cies in the flame (CH 4 , N 2 , 0 2 , H 2 6, CO, C0 2 , H 2 ). 
Species -specific, temperature-dependentquenching 
rates were taken from Tamura et al. [S]. The heat of 
formation of CH* was set at 66.3 kcal/mole above 
that of ground-state CH, based on the energy of the 
spontaneously emitted photon. 

OH* formation was modeled with a single reac- 
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Fig. 3. OH* mole fraction profiles determined bv measurement, GRI Mech 2,11, and an alternate hydrocarbon 
mechanism. 


Ill) CH + 0 2 <-> OH* + CO 

The rate constant for the reaction CH + 0 2 — > All 
Products was measured by Berman [16]. This mea- 
surement was used for k in by Marchese et al. [17] 
and will be used here as well. Uncertainty in this 
rate constant is likely to be large but is not readily 
estimated. Spontaneous emission and collisional 
quenching reactions were added to the mechanism 
in the same manner as previously described. The 
heat of formation of OH* was set at 93 kcal/mole 
above that of ground-state OH. Note that chemi- 
excitation reactions I — III all have other channels that 
form ground- state species, which are included in 
both kinetic schemes. 


Results and Discussion 

As with our previous investigations, GRI Mech 
2.1 1 continues to produce very good agreement with 
the overall flame length observed in the experi- 
ments, while significantly underpredicting the flame 
lift-off height. The alternate kinetic scheme is much 
more accurate in predicting lift-off height but ov- 
erpredicts the overall flame length. Note that the 
predicted lift-off height is related to the extinction 
strain rate obtained in the corresponding counter- 
flow diffusion flame. The extinction strain rate com- 
puted with GRI Mech 2.11 is nearly 20% higher 
than that obtained with the alternate mechanism, 
and thus the flame can anchor itself in a region of 
higher strain. 

The character of the CH distribution within the 
GRI-computed profile is in excellent agreement 


with the measurements; that is, the highest CH con- 
centration appears near the flame anchoring region, 
falling off to a nearlv constant level throughout the 
flame front, up to and including the tip of the flame 
front. The peak measured ground-state CH number 
density was 1.53 X 10 13 cur 3 . Given that the tem- 
perature was measured in this flame previously; we 
can easily calculate mole fractions from measured 
concentrations. Ilence, the measured peak mole 
fraction of CH is 4 X 10 f \ nearly 20% higher than 
the prediction of GRI Mech 2.11. 

Other than a better prediction of tbe lift-off 
height, the computed CH results are not as good 
with the alternate hydrocarbon kinetic mechanism. 
The location and distribution of CH is still within 
reasonable agreement, but the CH profile no longer 
appears to close at the flame tip, Additionally, the 
peak CH concentration is predicted to be 42% lower 
than the measured concentration. The dependence 
of the spatial characteristics and absolute concentra- 
tions on the choice of kinetic scheme reveals CH to 
be an important test of our rather well-validated 
flame model, as subtle differences are not seen be- 
tween the two kinetic schemes for the prediction of 
major species and temperature profiles. 

The CH* profile is thinner than that of ground- 
state CH when comparing measured CH to mea- 
sured CH*, as seen in Fig. 2. The difference in peak 
concentration of the experimental profiles of CH 
and CH* is of order 1000, similar to the ratio mea- 
sured by other researchers [IS]. A relationship simi- 
lar to that discussed with CH exists between the 
computed profiles of CH* and the measurements, 
as far as spatial distribution and overall characteris- 
tics. However, the measured peak mole fraction ot 
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TABLE 1 

Reactions added for CH* and OH* kinetics. Units for rate constants are centimeters, moles, and seconds, with the 
formulation k = AT B exp[ -E„/RT] - E a in cal niole^ 1 and R in cal mole 1 K 1 


Reaction 

A 

B 


Ref. 

C2H + O = CH* + CO 

LOSE + 13 

0.00 

0 

[15] 

C2H + 02 = CH* 4 - C02 

2.17E + 10 

0.00 

0 

[15] 

CH* -> CH 

1.85E 4- 06 

000 

0 

[8] 

CH* + N2 = CH + N2 

3.03E + 02 

3.40 

-381 

[8] 

CH* + 02 = CH -l- 02 

2.48E + 06 

2.14 

-1720 

[8] 

CH* + H20 = CH + H20 

5.30E 4 - 13 

0.00 

0 

[8] 

CH* + H2 = CH + H2 

1.47E 4- 14 

0.00 

1361 

[81 

CH* + C02 = CH + C02 

2.40E - 01 

4.30 

-1694 

[8] 

CH* + CO = CH + CO 

2.44E + 12 

0.50 

0 

[8] 

CH* + CH4 = CH + CH4 

1.73E + 13 

000 

167 

18] 

CH + 02 = OH* + CO 

3.25E 4- 13 

0.00 

0 

[16] 

OH* — » OH 

1.45E 4- 06 

0.00 

0 

[8] 

OH* + N2 = OH + N2 

LOSE 4 - 11 

0.50 

-1238 

[8] 

OH* + 02 = OH + 02 

2.10E 4- 12 

0.50 

-482 

18] 

OH* + H20 = OH + H20 

5.92E 4 - 12 

0.50 

-861 

18] 

OH* + H2 = OH + H2 

2.95E + 12 

0.50 

-444 

[8] 

OH* + C02 = OH + C02 

2.75E 4- 12 

0.50 

-968 

[8] 

OH* + CO = OH + CO 

3.23E 4- 12 

0,50 

-787 

[8] 

OH* + CH4 = OH + CH4 

3.36E + 12 

0.50 

-635 

[8] 


CH*, 2 X 10 " 9 , is a factor of 17 and 33 greater than 
the predictions of the alternate mechanism and GRI 
Mech 2.11, respectively. 

As with CH and CH*, the highest concentration 
of OH* is found at the anchoring region, as seen in 
Fig. 3. The measured OH* concentration profile ap- 
pears somewhat broader than predicted by either set 
of kinetics, and the measured peak mole fraction is 
within 15% of the peak computed with the alternate 
mechanism. 

Conclusions 

In this study, we extended the results of previous 
combined numerical and experimental investiga- 
tions of an axisym metric laminar diffusion flame in 
which difference Raman spectroscopy, laser-induced 
fluorescence, and a multidimensional flame model 
were used to generate profiles of the temperature 
and major and minor species. We discussed issues 
related to the computation and measurement of CH, 
CH*, and OH* in an unconfined laminar flame in 
which a cylindrical fuel stream is surrounded by a 
coflowing oxidizer jet. Experimentally, CH radical 
concentrations were measured with laser-induced 
fluorescence, whereas CH* and OH* concentrations 
were measured with flame emission. 

The results of this study indicate that GRI Mech 
2.11 does an excellent job of predicting peak CH 
concentration, considering the ppm concentration 
levels and narrow spatial extent of the CH radical’s 


profile, as well as the overall characteristics and 
shape of the CH profile. Other than a better predic- 
tion of flame lift-off height, the CH results are not 
nearly as good with our alternate hydrocarbon ki- 
netic scheme. 

As far as spatial distribution and overall character- 
istics are concerned, relationships similar to those 
observed between calculated and measured CH 
were observed for calculated and measured CH* 
and OH*, as well. Peak concentration levels for CH* 
were severely underpredicted with both kinetic 
schemes, while peak concentration levels of OH* 
agreed to within 15% (alternate mechanism) and a 
factor of 2 (GRI Mech 2.11) of the predicted peaks. 
Still, this indicates that although overall flame struc- 
ture is well understood and well characterized, the 
formation and destruction kinetics associated with 
excited-state species in flames requires further re- 
search. The reasonable uncertainties in the excita- 
tion reaction rate constants for CH* production in- 
dicate it is likely that new formation pathways, such 
as C 2 4- OH — > CH* 4- CO [11], need to be pos- 
tulated and investigated. 
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COMMENTS 


David R. Crosley, SRI International, USA. You predicted 
the CH concentration quite well, and of course you know 
CL even better. Thus, the correspondence between mea- 
sured and predicted OH* depends only on a single rate 
coefficient and suggests some adjustment is needed. Does 
the predicted OH*, which differs between the two mech- 
anisms, track directly the predicted CH? 

Authors Reply. The peak mole fractions of both CH and 
OH* were greater using GRI Mech 2.11 than with the 
alternate mechanism (see table). However, if one forms the 
ratio of peak CH to peak OH* mole fractions using the 
two different mechanisms, the ratios differ by nearly a fac- 
tor of 2. The location of the peak CH mole fraction in the 
GRI Mech calculation corresponds to a region with higher 
Oo mole fraction compared to the calculation using the 
alternate mechanism. In regions where the CH exists, the 
CO mole fractions and temperatures are nearly the same 
for both calculations. As a result, when the effects of all 
four species are taken into account, the ratio X <:n Xq/ 
Xoir Xco ls essentially the same for the two mechanisms. 


Peak Mole 

Fraction GRI Mech 2.11 

\ H 3.30E-06 

Xoh* 2.70E-08 

• 

Steve Hasko, BG Technology pic, UK. 

1 . Did the authors use die same rate data in dieir alternate 
reaction scheme as they used in their implementation 
of GRI Mech? Otherwise, what was that scheme opti- 
mized against? 

2. On the experimental and computed maps of tempera- 
ture, the tip of the inner cone of the computed result 
seems to be sharper than that of the measured flame, 
which appears more rounded. It has been our experi- 
ence that the strength of the tip can affect the lower 
edges and, thus, also the standoff height. Have the au- 
thors any thoughts about what might be happening at 
the tip, such as thermal diffusion of M atoms, that might 
cause this difference between the measured and com- 
puted results? 


Alternate 

Mechanism 

2.30E-06 

1.10E-08 



CH, CH*, OH* IN LAMINAR FLAME 


623 


Authors Reply. The two mechanisms did not use the 
same rate data. The alternate mechanism was compared 
against experimental measurements for temperature and 
species in a variety of flame configurations (Ref. [3] in die 
paper and references therein). 

The apparent sharpness at the tip of the inner cone seen 
in the computed temperature maps is due, in part, to the 
lower spatial resolution of the calculations at this down- 
stream location. We note that the calculations using the 
two different mechanisms display a similar degree of sharp- 
ness, despite having quite different lift-off heights. Your 
observation is an interesting one, but we have not had an 
opportunity to investigate this further. 


H. F. Calcote, Chernlori , Inc., LSA Ions in hydrocar- 
bon-oxygen flames are due to the che mi-ionization reac- 
tion: CH + O — » HCO + where the CH is usually con- 
sidered to be the excited state CH(A 2 A). Your work is thus 
of relevance toward quantitatively accounting for chemi- 
ions. The rate coefficient for the foregoing reaction with 
CH(A 2 d) has been reported as 4.8E + 14 [1] as measured 
in a flame by Cool and Tjossem [2J. How would the inclu- 
sion of this reaction in Table I affect your results? 

Table 1 includes three reactions with negative activation 
energies. Could you comment on these reactions in terms 
of Benson and Dobis recent discussion [3] of negative ac- 
tivation energy reactions? They argue that such reactions 
are an artifact of the experiment or they are due to a mul- 
tistep transition-state mechanism. 

REFERENCES 

1. Filakov, A. B., Prog. Energy Combust Sci. 23:399-528 
(1997). 

2. Cool, T. A. and Tjossem, P. J . H ., Chem. Phys. Lett. 
111:82-88 (1984). 


3. Benson, S. W. and Dobis, O., }. Phys. Chem. A 
1 02:5175-45 181 (1998). 

Author’s Reply. The rate constant for the che mi -ioniza- 
tion step is approximately a factor of 4 (for CO and H 2 ) to 
an order of magnitude (for 0 2 , N 2 , and C0 2 ) higher than 
the rate constants we used in our CH* collisional deacti- 
vation steps in the region of the flame where CH* is pres- 
ent. However, considering the magnitude of the O-atom 
concentration compared with the concentration of the col- 
lisional partners in this region (e.g., N 2 ), the chemi-ioni- 
zation step will have little impact on the consumption of 
CH*. 

Benson and Dobis discuss the apparent nature of neg- 
ative activation energies for certain bimolecular metathesis 
reactions. However, the negative activation energies 
quoted for three of the CH* collisional deactivation steps 
(as well as the OH* collisional deactivation steps) are nei- 
ther an experimental artifact nor due to a multistep tran- 
sition-state mechanism. The three parameter pseudo-Ar- 
rhenius forms used by Tamura et al. (our Ref. [8]), are 
taken from the compilation of Heinrich and Stuhl [1], who 
determined best fits of their experimentally determined 
quenching rate constants grouped with a variety of previ- 
ously quoted literature values in alternate temperature 
ranges. Indeed, Heinrich and Stuhl noted the difficulty in 
fitting the observed temperature dependencies with phys- 
ically meaningful rate parameters. While the sum of two 
Arrhenius forms was observed to be a good representation 
of the date, the authors found no indication that the 
quenching reactions occurred via two paths with distinct 
activation energies and thus opted for the more general 
three -parameter pseudo-Arrhenius form. 
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